The clustering of polarity reversals of the geomagnetic field 
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Often in nature the temporal distribution of inhomogeneous stochastic point processes can be 
modeled as a realization of renewal Poisson processes with a variable rate. Here we investigate one 
of the classical examples, namely the temporal distribution of polarity reversals of the geomagnetic 
field. In spite of the commonly used underlying hypothesis, we show that this process strongly 
departs from a Poisson statistics, the origin of this failure stemming from the presence of temporal 
clustering. We find that a Levy statistics is able to reproduce paleomagnetic data, thus suggesting 
the presence of long-range correlations in the underlying dynamo process. 

PACS numbers: 91.25.Mf; 91.25.-r; 02.50.-r; 91.25.Cw 



Local paleomagnetic measurements of the geomagnetic 
field 0&0 reveal a sequence of sudden and occasional 
global polarity reversals in the last 160 million years. The 
magnetic dipole typically reverses its direction in a few 
10 3 years, while the time intervals between successive re- 
versals range from 10 4 up to 10 7 years IE 0- Even 
though polarity reversals can be ascribed to the Earth's 
magnetic dynamo IE IE j details of the mechanism 
remain poorly understood. The debate on the origin of 
reversals, the modeling of the trigger (external or inter- 
nal to Earth) and the generation of longer variations in 
the average reversal rate is still open (cfr. e.g. Ref.s 
13 )• A look at the sequence of reversals, for exam- 
ple the most used CK95 database 0, reproduced in fig. 
^ shows that polarity reversals seem to be the result of 
a non-periodic chaotic (or stochastic) process. Actually 
non regular reversals can be observed in the framework of 
purely deterministic toy models that mimic the dynamics 
of the dynamo effect with only few modes [E HE HjJ > m 
the framework of noise-induced switchings between two 
metastable states |0, ^E 13 or in mean-field dynamo 
with a noise-perturbed a profile El- In principle 
geodynamo is described by 3D global Magnetohydrody- 
namics (MHD) that self-consistently solve for the fluid 
flow, thermodynamics and magnetic fields with all non- 
linear feedbacks (for a review see Ref. and references 
therein, and the results of the recent 15.2 TFlops simu- 
lation of geodynamo on the Earth Simulator |l(j)- Al- 
though some numerical codes have simulated short series 
of spontaneous reversals, none have been run so far at 
high enough resolution to be confident that the critical 
dynamics is being captured in the simulation. 

In spite of the paucity of data sets, it is commonly as- 
sumed that the phenomenon of reversals stems from an 
underlying Poisson process. This conjecture is based on 
the fact that the polarity persistence times At, defined 
as the time intervals between two consecutive reversals 
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FIG. 1: Upper panel: probability density function P(At) of 
persistence times At (statistical errors are shown as vertical 
bars). Lower panel: polarity of the earth's magnetic field 
(from today). The solid bar are the normal (present) polarity. 
We used the CK95 dataset. 



At = tj+i — ti, seem to be exponentially distributed 
EE EEC!; namely P(At) ~ Aexp(-AAi) where A 
represents the rate of occurrence of reversals. Different 
analyses (cfr. e.g Ref. ^3) state that the frequency 
distribution of intervals between Cenozoic geomagnetic 
reversals approximate a power law for large At. Even a 
rough look at the probability distribution function P(At) 
(cfr. fig^) shows that the situation is not clear, mainly 
in presence of a poor data set with high statistical errors. 
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FIG. 2: The probability densities f(H) of the stochastic vari- 
able H (left panel) and the cumulative probability P(H > h 
(right panel), are presented for all datasets. Theoretical prob- 
abilities observed under a Poisson statistics are also shown. 



of a Local Poisson Process is correct or not. In other 
word, according to Constable 0] it can be conjectured 
that (hypothesis H ), even if it cannot be decide wheter 
globally the reversals stem from a Poisson process, an un- 
derlying time-varying Poisson process could be assumed 
to originate the geomagnetic reversals. Since the rever- 
sals rate A(t) is not known, we have to build up a test 
which is independent on the rate A. This can be done 
through a measure used previously in the framework of 
solar flares 0, HH< We introduce a statistical quan- 
tity that is nothing but the suitably normalized local 
time interval between reversals. Let us consider the time 
sequence of reversals as a point-like process and sup- 
pose that each reversal occurs at the discrete time U. 
Let Sti = min{ii + i — tf,ti — f,_i} and let Tj be either 
n = ti-i - tj-2 (if Sti = U - U-i) or Ti = t i+2 - U+i 
(if Sti = ti+i — U), so that Sti an d are the two per- 
sistence times following or preceeding a given reversal at 
tj. If the local Poisson hypothesis Ho holds, both 6U and 
Ti are independently distributed according to an expo- 
nential probability density: p(St) = 2Aj exp(— 2\i$t) and 
P(t) = Ai exp(— A»t) with local rate Aj. Then, under the 
hypothesis Hq, the stochastic variable H, defined as 



Recently Constable [Tg raised two main statistical 
features to the attention of the scientific community, 
namely: i) the distribution of events shows a paucity 
of short intervals; ii) the rate of occurrence of events 
is time-dependent A = A(f). The author investigated 
the temporal distribution of A(t), showing that a Pois- 
son model with a monotonic rate, either increasing or 
decreasing, is not a good model for the reversal process. 
On the contrary, reversals could be perhaps modeled as 
a renewal Poisson process with a rate that must change 
sign at some interval before 158 My In any case, 

modeling A(t) over the entire time interval is not an easy 
task (cfr. also Ref.s [lH l2*l|L Moreover, in this situa- 
tion, namely when the rate A depends on time, the dis- 
tribution of persistence times remains without physical 
meaning and cannot be used to determine the Poisson 
character of events 20] . 

The problem is perhaps more general than what we 
present here, because abrupt flow reversals have been 
found also in the large-scale circulation during turbulent 
Rayleigh-Benard convection j^, or in the wind direc- 
tion in atmosphere [23|. The conjecture that reversals 
are Poisson events is made in all cases. Here, starting 
from the above experimental evidences, and using a sim- 
ple statistical test on some databases, we will investigate 
whether a conjecture based on the occurrence of a Pois- 
son process for reversals is correct or not. We will show 
that this is not the case, and that geomagnetic rever- 
sals are clustered in time, a result which stems from the 
presence of memory in the process generating polarity 
reversals. 

Under the experimental evidence that the rate of re- 
versals is not constant, we can test, as a zero-th order 
hypothesis, wheter an approach based on the occurrence 



H(Sti 



2SU 



2St % 



(1) 



is uniformly distributed in [0; 1] and has cumulative dis- 
tribution P(H > h) = 1 — h, independent on Xi (2j|. 
In a process where r^s are systematically smaller than 
2<5tiS, clusters are present and the average value of H is 
greater than 1/2. On the contrary when the underlying 
process is characterized by voids, the average value of H 
is less than 1/2. From time series, it is easy to calculate 
the probability P(H > h) and the probability density 
function f(H). 

We apply the above test to four different sequences 
of geomagnetic polarity reversals, namely to the already 
mentioned CK95, to H68, HA97 and KG86 Hj. We cal- 
culate the probability density function f(H) reported in 
hgure (left panel). As can be seen, a significant de- 
viation from the uniform distribution is evident in all 
datasets, the departure of polarity reversals from local 
Poisson statistics stemming from the presence of clusters. 
Then a clear deviation of the observed cumulative proba- 
bility P(H > h) from a linear law (cfr. fig. fright panel), 
expected under Hq, is obtained. A Kolmogorov-Smirnov 
(KS) test applied to the cumulative distributions con- 
firms that the assumed hypothesis Hq is not reliable and 
must be rejected (the significance level of the KS test 
being smaller than 0.5% for all datasets). 

To further characterize the origin of the departure of 
polarity reversals from a local Poisson process, we try to 
describe the statistics of persistence times with a Levy 
process. The Levy functions are stable, and can be ob- 
tained from the Central Limit Theorem by relaxing the 
hypothesis of finite variance j2^. To avoid problems 
arising from the infinite variance, which is unlikely in 
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FIG. 3: The cumulative probability P(At > T) as a function 
of T (symbols) obtained from the CK95 database (vertical 
bars represent statistical errors). The full line represents the 
fit obtained with the truncated Levy function ©. The best- 
fit parameters are v = 1.09 ± 0.1 and 6 = 0.15 ± 0.04. 

real physical processes, a Truncated Levy Flight distri- 
bution (TLF) has been proposed |2!j introducing a cutoff 
in the standard Levy. The process can be generated by a 
random variable x according to the distribution function 
f{x) ~ e 61 1 ^ I | x I" 1 "", where 9 represents the cutoff rate 
and v is the characteristic exponent, ft can be shown |3(| 
that the limit distribution P(z) of the sum z of random 
variables x can be computed, for v ^ 1, by the following 
integral 

f°° f 9 U 

P(z) — CI dk cos(zfc) exp < — — 

K ' Jo P \cos(7nV2) 



2n(z 2 + 2 y/ 2 cos[v arctan(| z \ /<?)] 
^r(^) sin(7rz/) 




where C is a normalization factor and T(v) is the usual 
gamma function. For v > 2 we recover a normal random 
process, the result of the integration being a Gaussian 
PDF, while smaller values of v represent increasing de- 
viation from Gaussianity. We thus fitted the measured 
cumulative PDF of persistences P(At > T), computed as 
described above and reported in figure [21 with equation 
(J2J. The equation has been numerically integrated us- 
ing a standard minimum-^ 2 procedure. We obtained the 
best-fit parameters v = 1.09 ± 0.1, and 9 = 0.15 ± 0.04 
with a reduced \ 2 — 0.5. 

The possibility of reproducing with a Levy function 
the cumulative distribution of persistence times indicates 
both that the process underlying the polarity reversals 
is statistically self-similar in time, and that a certain 
amount of memory, due to long-range correlations, is 
present in the process. From a physical point of view 
we could expect that the dynamics of the fluid earth core 
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FIG. 4: The probability densities f(H) of the stochastic vari- 
able H (left panel) and the cumulative probability P(H > h 
(right panel), are presented for both models A and B. The- 
oretical probabilities observed under a Poisson statistics are 
also shown. 



is affected by its hystory, thus generating correlations 
among reversals. This feature has recently be observed 
for the solar dynamo jUl- Of course a lot of different 
random processes can reproduce the departure from a 
Poisson statistics. Only the next generation of simula- 
tors will be able to produce datasets with enough rever- 
sals to allow us toinvestigate in detail the occurrence of 
long-range correlations in 3D global MHD simulations of 
geodynamo. In this framework, statistical analysis on 
real data, among other, play the key role of discriminat- 
ing among different random processes that can reproduce 
the departure from poisson statistics, thus increasing our 
knowledge of the geodynamo process. Dynamical models 
can help us because they describe, with only few physi- 
cal ingredients, some gross features of the enormous com- 
plexity of the geodynamo process. In this perspective it 
is useful to compare the statistics of reversals observed 
in toy models with statistics obtained on real datasets. 

As an example, we investigate the sequence of random 
reversals generated by a standard numerical analysis of 
the two-disk chaotic geodynamo model 0- The model, 
known as "Rikitake dynamo" , is described by the follow- 
ing ordinary differential equations: dx/dt = —x + yz, 
dy/dt = —y + x(z — 15/4), and dz/dt = 1 — xy (in the 
following model A). The cumulative distribution of H, 
obtained by the times of reversals of the variable x(t), has 
been reported in figure 01 It is evident that the chaotic 
dynamics of the model is responsible for the presence of 
correlations among the reversals. As a further example 
we investigate a dynamical model based on stochas- 
tic excitation of the axisymmetric component of the mag- 
netic field (model B), that reads: dx/dt = (1 — x 2 )x + 
Vnx + Vi 2 y + Vi 3 z, dy/dt = -ay-cz + V 2 ix+V22y+V 23 z, 
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and dz/dt = cy — az + V31X + V^y + V33Z. The Vij(t) 
are 9 independent random functions of time, with zero 
mean and equal r.m.s. magnitude, that are renewed af- 
ter a time r c . We use the standard parameters a = 2, 
c = 5, r c = 0.01, and D = (V*) r c , where D = 0.4 [l4|. 
The cumulative distribution of H (fig. 0J, calculated by 
using the times of the random reversals of the variable 
x(t), has been reported in figure^ Even in this case a 
departure from a Poisson distribution is observed, due to 
the correlations introduced through r c . 

Even if the chaotic dynamics within the Rikitake 
dynamo roughly reproduces the presence of clustering 
among reversals, a look at fig.s0]reveals that in model A 
the departure from the local Poisson distribution seems 
to be affected by the presence of sudden jumps in the 
cumulative distribution of H. Say the reversal time se- 
ries x(t) presents few recurrent persistence times of equal 
extent. On the other hand the model B shows a differ- 
ent departure from Poisson statistics, that is mainly due 
to both very low values and, even if in minor extent, to 
very high values of H . The time behaviour of reversals in 
the model B is then hardly compatible with the cluster- 
ing process evidenced through our analysis. Further toy 



models based on completely random external triggers of 
reversals (e.g. stochastic resonance) cannot describe the 
geodynamo process. 

As a conclusion, we investigated the statistics of persis- 
tence times between geomagnetic reversals. We applied a 
statistical test, showing that geomagnetic reversals stem 
from an underlying process that is far from being locally 
Poissonian, as currently conjectured |Tsj . A Levy func- 
tion is able to nicely fit the probability distribution of 
persistence times. Although not investigated up to now 
in a geophysical framework, our results can be interpreted 
as a strong evidence for the presence of correlations be- 
tween reversal events. These correlations arise from some 
de gree of memory in the underlying geodynamo process 
[71 l27j that gives rise to clustering of reversals. 
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